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ABSTRACT 

Context. High redshift quasars can be used to deduce the distribution of dark energy in the Universe, as a complementary tool to SN 
la. 

Aims. The method is based on determination of the size of the Broad Line Region from the emission line delay, determination of the 
absolute monochromatic luminosity either from the observed statistical relation or from a model of the formation of the Broad Line 
Region, and determination of the observed monochromatic flux from photometry. This allows to obtain the luminosity distance to a 
quasar independently from its redshift. The accuracy of the measurements is however, a key issue. 

Methods. We model the expected accuracy of the measurements by creating artificial quasar monochromatic lightcurves and responses 
from the Broad Line Region under various assumptions about the variability of a quasar. Broad Line Region extension, distribution 
of the measurements in time, accuracy of the measurements and the intrinsic line variability. 

Results. We show that the five year monitoring based on Mg 11 line should give the accuracy of 0.06 - 0.32 magnitude in the distance 
modulus which allows to put interesting constraints on the cosmological models. Monitoring of higher redshift quasars based on CIV 
lines is problematic due to much higher level of the intrinsic variability of CIV in comparison with Mg II. 

Key words, accretion, accretion disks - black hole physics 



1. Introduction 



The key issue of the present day cosmology is the determina- 
tion of the properties of the dark energy distribution (for a re- 
view, see Frieman, Turner & Huterer 2008). Watson et al. (2010) 
proposed that Active Galactic Nuclei (AGN) can also be used 
as probes of this distribution, complementary to the study done 
with other methods. The basic idea is the following. The ab- 
solute size of the Broad Line Region (BLR) in AGN can be 
measured from the time delay between the variable continuum 
emission from the nucleus and the emission lines in the opti- 
cal/UV band (e.g. Blandford & McKee 1982; Peterson 1993; 
Sergeev et al. 1994; Wandel et al. 1999). This reverberation 
technique based on the Hy? line was succesfully used for low 
redshifts objects (up to z = 0.292; Kaspi et al. 2000, Bentz et 
al. 2009). The size of the BLR is related to the absolute value 
of the source monochromatic luminosity. This was found ob- 
servationally (Kaspi et al. 2000; Peterson et al. 2004; Bentz et 
al. 2009), but the relation was later theoretically explained by 
Czerny & Hryniewicz (201 1) as corresponding to the formation 
of the inner radius of BLR at fixed temperature of 1000 K in 
all sources, independently from the black hole mass and accre- 
tion rate. This results points toward the important role of dust 
in the accretion disk atmosphere in formation of BLR, connect- 
ing it with the dust sublimation temperature. Having the absolute 
value of the monochromatic luminosity, and measuring the ob- 
served monochromatic luminosity, we can determine the lumi- 
nosity distance to the source. Summarizing, the measurement of 
the observed monochromatic flux and of the time delay between 
continuum and an BLR emission line for a number of sources 



allows to construct the diagram of the luminosity distance as a 
function of redshift, and to test the cosmological models. 

The use of an independent tracer of the dark energy is impor- 
tant. Various methods are being used or will be used to trace the 
dark energy, but each of them can have certain bias. Supernovae 
luminosity may be subject to systematic cosmological evolution, 
due to the evolution of the chemical content of the Universe. 
Quasars show roughly solar or even super-solar metallicity even 
to redshifts of a few (e.g. Dietrich et al. 2003a,b; Simon & 
Hamann 2010; De Rosa et al. 2011; Gafl et al. 2011), proba- 
bly due to the correlated evolution of starburst and quasar ac- 
tivity. Therefore, the relation between the size of the BLR and 
monochromatic luminosity, reflecting dust properties, should not 
contain strong evolutionary trends with the redshift. Other meth- 
ods, like galaxy clusters, weak gravitational lensing and galaxy 
angular clustering (baryon acoustic oscillations) can also con- 
tain some bias since the clusters are not spherically symmetric 
and relaxed, and galaxy clustering may be the subject of magni- 
fication bias. If numerous methods are used, the results will be 
more reliable. 

We have recently started an observational campaign with the 
South African Large Telescope (SALT) to collect data for re- 
verberation studies of a sample of 7 objects. We are planning to 
use the low ionization line of Mg II 2798A, which is suitable 
for studies of objects in the redshift range from 0.9 to 1.4. So 
far only a few authors claimed the measurement of the delay in 
distant quasars using Mg II or CIV line (Wang et al. 1997 for 
3C 390.3, z = 0.0561, using CIV line; Chelouche et al. 2012 
for MACHO 13.6805.324, z = 1.72, from photometric method). 
Woo (2008) showed that using Mg II line is a viable possibility 
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but the two year monitoring program was of course too short to 
attempt any delay measurement. Reverberation studies of even 
higher redshift quasars, based on the CIV 1549 A line, were re- 
cently done with the HEX telescope by Kaspi et al. (2007), but 
the authors generally were not successful in the measurement 
of the time delays, with a tentative delay measurement of 595 
days for one object (S5 0836+71) out of six. Therefore, in this 
paper we perform a detailed analysis of the requirements for a 
successful monitoring program. 



2. Method 

AGN are highly variable in all energy bands (for a review, see 
Ulrich et al. 1997), including optical band (e.g. Matthews & 
Sandage 1963; Hawkins 1983; Barbieri et al. 1988; Cristiani et 
al. 1996; Vanden Berk et al. 2004; Bauer et al. 2009; MacLoad 
et al., 2012). Although best studies were done in the X-ray band, 
there were also several attempts to describe quantitatively the 
character of the optical variability. The power density spectra 
(PDS) in the optical band for NGC 4151 were constructed by 
Terezbizh, Terebizh & Biryukov (1989) and by Czerny et al. 
(2003) in the long time-scale, and by Edelson et al. (1996) in 
the short time-scale; for NGC 5548 by Czerny et al. (1999); 

Therefore our idea here is to create simulated datasets of a 
fiducial high redshift quasar to estimate the role of various pa- 
rameters of an observational campaign (e.g., number and sepa- 
rations of observations, duration of the campaign) and expected 
quality of data, in the accuracy of the determined reverberation 
time delay. This should help in developing the optimum strategy 
of such reverberation studies campaign. 

We use the Timmer & Konig (1995) algorithm for generating 
an artificial continuum lightcurve of a quasar based on assumed 
shape of the power density spectrum (PDS). We then expo- 
nentiate the resulting lightcurve, following Uttley et al. (2005), 
since this reproduces the log-normal distribution characteristic 
for high energy lightcurves of accreting sources. We assume that 
the PDS has a shape of a doubly-broken power law, although the 
long term variability is highly uncertain, but this allows us to test 
whether there are any effects of the power leak toward shorter 
frequencies. We fix the longest frequency slope at zero, while 
the two remaining slopes, si and S2, and the break frequencies, 
/i and /2, are in general the free parameters of the simulations. 
We expect that they are mostly affected by the value of the black 
hole mass of the monitored quasar, since in the X-ray band the 
black hole mass is the key parameter of the PDS (Hayashida 
et al. 1998; Czemy et al. 2001; Papadakis 2004, Zhou et al. 
2010; Ponti et al. 2012), although the efi'ect of the Eddington 
ratio and/or spectral state is an issue (e.g., McHardy et al. 2006; 
Nikolajuk et al. 2009). In most cases we fix these parameters at: 
si = 1.2, S2 = 2.5, /i = 0.012yr-i and/2 = 0.18yr ^ These val- 
ues were estimated from the optical PDS of NGC 5548 (the in- 
termediate slope, si and the break frequencies, /i and f2), while 
the high frequency slope, S2 was taken from the Kepler quasars 
monitoring campaign, after correcting the slope for the subtrac- 
tion of the Unear trend (Mushotzky et al. 2011). The break fre- 
quencies were then scaled to a higher expected black hole mass, 
and for the redshift. The overall dimensionless amplitude of the 
PDF is also a model parameter, but again in most simulations it 
was fixed at 0.1 and 0.3, which correspond to the dimensionless 
dispersion in the infinite photometric lightcurve (i.e., of duration 
above 1 000 yr) without any measurement error of 0.32 and 1.61, 
correspondingly. The created lightcurves have the total standard 
length of 30 years and are densly spaced (1 day). 



Next we model the response of the BLR (meant to be repre- 
sented by the Mg II 2798A line) to the variabiUty of the contin- 
uum. The line lightcurve is created by convolving the continuum 
lightcurve with a gaussian kernel. The center of the gaussian rep- 
resents the Une delay after the continuum, which is a free param- 
eter of the model, and is of the order of a few hundred days. The 
gaussian width models the smearing of the response due to the 
finite extension of the BLR. Most of the reprocessing is done by 
the BLR part located on the far side of the black hole since the 
clouds on the side to the observer are dark and so they do not 
emit efficiently. Therefore, in most computations we assume the 
gaussian width of the order of 10 percent of the time delay, since 
such parameters recover well the measured accuracy of the time 
delay for NGC 5548. 

Obviously, the photometric and spectroscopic monitoring 
with SALT telescope cannot be done with 1 day spacing, so in 
the next step we select observational points from the two light 
curves, assuming a realistic observational pattern and simulat- 
ing also the measurement errors. In most cases we assume that 
the object cannot be observed for about 3 months every year (it 
is then out of the observational window of SALT), and the gap is 
the same for photometric and spectroscopic lightcurve. For pho- 
tometric monitoring we assume that the object will be observed 
roughly every two weeks, with random dispersion of 3 days. For 
spectroscopic measurements, we assume 5 SALT observations 
every year, with random dispersion of 14 days. We can of course 
vary these parameters as well. 

Finally, we have to account for possible uncertainties of both 
the spectroscopic and photometric measurements. Therefore, the 
value of the Ughtcurve at each observational point is replaced 
with a random value obtained from a gaussian distribution with 
the mean equal to this value and the dispersion corresponding 
to the assumed 1 sigma error. The errors for the line measure- 
ment are generally assumed to be larger than the errors of the 
photometric measurements. 

Having the two curves, we now measure the time delay 
of the spectroscopic lightcurve with respect to the photometric 
one, for example using the interpolated cross-correlation func- 
tion (ICCF) as usually done for the real data (e.g. Gaskell & 
Peterson 1987). We interpolate the spectroscopic measurements 
to the dates of the photometric measurements, and we use the 
maximum of the ICCM to determine the delay. In order to deter- 
mine the expected accuracy of the measurement, for each set of 
the all other parameters we analyze 1000 random realizations of 
the process, and we create the histogram of the measured delays. 
The comparison of this histogram with the assumed delay char- 
acterizes the expected quality of the observational campaign. In 
particular, we use the median of the histogram to give the best 
value of the delay, and for the dispersion we give the dispersion 
around the known expected value. In addition we give the spread, 
defined as a half range between the 1-0.158 and 0.158 quantiles 
of the histogram. For the Gaussian distribution it would corre- 
spond to Icr error. We also experiment with alternative methods 
to measure time delay. 

The whole procedure should mimic well the actual reverber- 
ation measurements performed for numerous less distant sources 
(e.g., Kaspi et al. 2000; Peterson et al. 2002; Bentz et al. 2009; 
Denney et al. 2009; Walsh et al. 2009; Denney et al. 2010; Barth 
et al. 2011; Stalin et al. 2011; Dietrich et al. 2012; Grier et al. 
2012). 

With our Monte Carlo approach we are then able to study 

which elements of the campaign are the most important ones 
and where there are any possibihties for an improvement. These 
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we hope will be of significant help in planning a c.a. 5 years long 
observational project. 

3. Results 

We model the expected accuracy of the measurement of the 
emission line delay with respect to the continuum for a single 
object (high redshift quasar) monitoring. The problem is im- 
portant because the variability of AGN are of a red noise type, 
with a broad range of timescales and a strong contributions from 
the timescales much longer than the observed lightcurves. The 
Monte Carlo approach is effective in assessing the likely results 
as well as possible caveats since in the real monitoring we have 
to work just with one realization of the process. 

The basic model (hereafter model A) is adjusted to our al- 
ready started observational campaign. The feasibility study has 
been done for a quasar Q 2113-4538 (z - 0.946). The high red- 
shift of the source ensures that the Mg II 2798A line is within 
the spectroscopy frame. The quasar comes from the Large Bright 
Quasar Sample (LBQS; e.g., Hewett et al. 1995), the black hole 
mass in this object was estimated to be 1 .5 xIO'^Mq (Vestergaardt 
& Osmer 2009), and the expected time delay of the Mg II line 
(from the Kong et al. 2006 formula) is 660 days, if the addi- 
tional increase due to the redshift is included. The spectroscopy 
is being done with the SALT and the source is visible for about 
9 months. We plan to perform 5 observations per year and two 
observations were already performed. The photometry will be 
performed more frequently, with smaller telescopes, approxi- 
mately every 2 weeks. For the first year, the collaboration with 
the OGLE team has been secured. 

We have already analyzed the first two spectra obtained with 
the SALT telescope, one in Dark Moon conditions, and one in 
the Bright Moon conditions, and the statistical error in the Mg II 
2798Aline equivalent width (EW) of the two spectra were 0.3%. 
However, possible uncertainty in determination of the level of 
the Fe II pseudo-continuum can contribute to the estimate of the 
Mg II line EW, and some other systematic errors are also likely 
to be present. Therefore we finally adopt 10% as the actual mea- 
surement accuracy of the line intensity in the basic model. The 
photometric accuracy of continuum measurement in the basic 
model is set to 5%. We give all the parameter values in Tables[T| 
and m 

Model A represents our reference setup of the campaign, 
while other models (from B to E) illustrate the dependence of 
the final results on the details of the campaign. We repeat the 
computations for two values of the amplitude of the continuum 
variability: a most likely one (first case) and a larger amplitude 
(second case). 

First we test the ICCF method to determine the time delay in 
the simulated data, since this method is most often used in real 
data analysis. 

Running 1000 realizations for model A we properly recover 
the assumed time delay of 660 days, but the median value has 
a noticeable error: 605 + 285 days for low variability case and 
605 + 268 for the other case. Such a measurement is not quite 
satisfactory, although if the monitoring is done for a number of 
objects, the statistical error in determination of the dark energy 
distribution will be reduced. However, there are also possibili- 
ties of reducing the error of the delay measurement for a single 
object with a more careful approach to the monitoring. 

First, we test the importance of the accuracy of measure- 
ments. Model B represents the situation when the individual 
measurement errors are significantly reduced, with the remain- 
ing organization of the campaign schedule unchanged. This, 



however, does not reduce significantly the uncertainty of the fi- 
nal result. Denser photometric or spectroscopic monitoring does 
not reduce the dispersion around the delay, either (see Table |2]i, 
if the other parameters remain unchanged. Only extending the 
monitoring to much longer timescales, of order of 20 years, re- 
sults in significant improvement. This is due to the strong red 
noise type variability of the source. For the assumed parame- 
ters, model A shows the dispersion in the photometric points of 
0.16 and 0.57, and the dispersion in the spectroscopic points of 
0.18 and 0.52 for the low and high variability level, respectively. 
Still, long timescale trends dominate the source variability and 
this leads to a significant fraction of the lightcurves with time 
delay measured very badly. We analyze this problem in detail 
below. 

3.1. ICCF problem 

The distribution of the delays obtained for 1000 realizations of 
the model A, low variability level is shown in Fig.[T] The distri- 
bution contains significant tail of measured time delays which 
are shorter than the assumed value. The presence of the tail 
skews the distribution (the mean delay is 549 days, while the 
median is 605), which enhances the expected error of the time 
delay measurement. If we simply neglect the solutions with the 
delay smaller than 500, then the resulting delay is 692 + 165; if in 
addition we skip the delays longer than 820 days, the estimated 
delay is 626 + 83, i.e. the error is much lower than previously, 
but we then skip 32.8% of the cases. For the higher variability 
level the longer delay tail is less important, but the shorter delay 
tail is very important. 

Thus with ICCF method, the tails of the distribution strongly 
affect the results. In the real situation we will have just one 
lightcurve for a single object, so there is a considerable prob- 
ability that the measurement for a single lightcurve will be com- 
pletely wrong. We therefore should characterize the probability 
of the proper delay recovery much better For a simple gaussian 
distribution, one parameter describes both spread of the distri- 
bution core and extent of its tails. The former relates to preci- 
sion of estimates, the latter describes reliability i.e. probability 
of a gross error. Delay statistic^ based on our simulations have 
strongly non-gaussian distribution, hence at least two parameters 
are needed to characterize their precision and reliability. For the 
latter characteristics we selected standard deviation (s.d.), for the 
former half spread of symmetric quantiles limiting central 68.3% 
of the distribution. 

777 

In order to illustrate better the source of such large errors we 
plot two examples of the lightcurve realization for model A, and 
the corresponding ICCF plots (see Fig.|2]i. The problem is even 
more severe if the overall variability amplitude is assumed to be 
lower although in this case the large outbursts happen less fre- 
quently. The ICCF still frequently returns spurious results when 
the lightcurves are dominated by a single monotonic event (flare) 
in the observed window. In addition, lower amplitude makes the 
observational errors relatively more important, so the overall ac- 
curacy of the delay determination is lower. 

This suggests a considerable space, and actual need, for im- 
provement. Therefore, we consider various alternatives to the 
ICCF method. To evaluate performance of the delay statistics 
we employ two criteria: its bias and the error Inspection of his- 
tograms obtained from our simulation reveals asymmetric distri- 
bution of the values of the statistics. Hence we prefer median 

' here statistics means a function of data, used for estimates 
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Table 1. Parameters of the exemplary models and the accuracy of the delay recovery for variability amplitude 0.1, assumed delay 
660 days 



Model 
Parameter 


A 


B 


C 




D 


E 


Photometric accuracy 


5 % 


1% 


5 % 




5 % 


5 % 


Spectroscopic accuracy 


10% 


1% 


10% 




10% 


10% 


Width of the BLR response 


10% 


10% 


10% 




10% 


10% 


No of photometric points/year 


18 


18 


18 




100 


18 


No of spectroscopic points/year 


5 


5 


5 




5 


18 


Dispersion of the photometry date 


3 days 


3 days 


3 days 




1 day 


3 days 


Dispersion of the spectroscopy date 


14 days 


14 days 


14 days 


14 days 


3 days 


Observational gap/year 


3 months 


3 months 


3 months 


3 months 


3 months 


Campaign duration 


5 years 


5 years 


20 years 


5 


years 


5 years 


ICCF 














Obtained delay 


605 ± 285 


605 ±251 


660 ± 94 


605 


±279 


605 ± 246 


Spread 


220 


192 


55 




220 


165 


x' 

Obtained delay 


671 ± 113 


652 ± 20 


673 ±31 


671 


± 111 


671 ±76 


Spread 


84 


9 


40 




93 


47 


ZDCF 

Obtained delay 
Spread 


622 ± 200 
141 


666 ± 38 
20 


661 ±405 
62 


604 


±200 
137 


641 ± 206 
119 



Table 2. Parameters of the exemplary models and the accuracy of the delay recovery for variability amplitude 0.3, assumed delay 
660 days 



Model 


A 


B 


C 


D 


E 


Parameter 












Photometric accuracy 


5 % 


1% 


5 % 


5 % 


5 % 


Spectroscopic accuracy 


10% 


1% 


10% 


10% 


10% 


Width of the BLR response 


10% 


10% 


10% 


10% 


10% 


No of photometric points/year 


18 


18 


18 


100 


18 


No of spectroscopic points/year 


5 


5 


5 


5 


18 


Dispersion of the photometry date 


3 days 


3 days 


3 days 


1 day 


3 days 


Dispersion of the spectroscopy date 


14 days 


14 days 


14 days 


14 days 


3 days 


Observational gap/year 


3 months 


3 months 


3 months 


3 months 


3 months 


Campaign duration 


5 years 


5 years 


20 years 


5 years 


5 years 


ICCF 












Obtained delay 


605 ± 268 


605 ± 260 


660 ± 126 


605 ± 261 


660 ± 247 


Spread 


220 


220 


28 


220 


165 


x' 

Obtained delay 


671 ±51 


652 ± 17 


673 ± 14 


652 ± 48 


652 ± 29 


Spread 


28 


9 


13 


28 


28 


ZDCF 












Obtained delay 


662 ± 106 


666 ± 39 


662 ± 30 


660 ± 104 


667 ± 112 


Spread 


41 


18 


26 


38 


38 



rather then average value for estimation of bias. For an arbi- 
trary distribution the former yields equal probability of under- 
estimated and overestimated result. For estimation of error we 
employ standard deviation with respect to the assumed delay 
value of 660 days. As it depends on quadratic norm of residu- 
als it is rather sensitive to gross errors. As we aim to minimize 
the gross errors, the standard deviation seems the right choice. 
However, to better characterize our resulting distribution we ad- 
ditionally list the half spread of their 1-0.158 and 0.158 quan- 
tiles. For the gaussian distribution it would coiTespond to half 
of the +lcr range, i.e. to Icr. A value of the spread small com- 
pared to the standard deviation indicates a distribution with the 
narrow core and broad wings. Note that generally for an asym- 
metric distribution center of the spread range does not coincide 
with the median or average. The corresponding values for all the 
setups are given in Table [T] and Table |2l 



3AA.x^ method 

In this method we shift the spectroscopic curve and interpolate 
to the photometric time measurements as before, but, instead of 
computing ccf, we minimize the function 



X 



1 ^ 
N 



(Xj - Ay if 



2' 



(1) 



where the constant A adjusting the normalization of the two 
curves can be obtained analytically from the condition of mini- 
mizing the function with respect to A: 



A = 



25 xyT, 



1/2 



(2) 
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0.008 



0.006 



0.004 



0.002 



0.000 



0.008 



0.006 



0.004 



0.002 



0.000 



0.008 



0.006 



0.004 



0.002 



0.000 



0.008 



0.006 



0.004 



0.002 



0.000 




ZD.CF.centrQid. 



ICCF. 







200 



400 600 
Lag [d] 



800 



Fig. 1. The normalized histogram of the time delays in 1000 re- 
alizations of the model A, low variability, for various methods: 
(upper panel), ZDCF and ICCF method (lower panel). For 
ZDCF we show results both for maximum and gaussian centroid 
method of delay determination, for other methods we give the 
maximum method. Altough the average value describes well the 
true delay insterted into the model, there is a considerable tail of 
much shorter or much longer values which contribute consider- 
ably to the overall dispersion. 



4 - lightcurve No, 1 



lightcurve No. 15 



12 3 

years 



3 3 

years 



500 1000 

delay [days] 



500 1000 
delay [days] 



Fig. 2. Two examples of the artificial photometric and spectro- 
scopic lightcurves and the corresponding plots of the ICCF, for 
model A. 



with the coefficients given by 



>xy3 - 



>x3y - 



i=l 
N 

i=l 
N 

^XiyiSxf. 



(3) 



The method works properly only in the case when the 
searched time delays are smaller than a half of the campaign 
duration. For longer delays the method picks up the longest val- 
ues of the time delay since this leaves only relatively short data 
sequences, where the variability is low due to the red noise char- 
acter of the curves, and then the two curves can be fitted very 
well independently of their actual shape. 

3.2. ZDCF method 

We also analyzed the delay between the simulated photometric 
and spectroscopic curves using the discrete correlation function 
(DCF). The idea is rather simple: all possible pairs are formed 
from two observation sets and then they are binned according to 
their time lag. The hitch is in the details of statistical evaluation 
of the bin mean and variance. Here we rely on the ZDCF flavor 
of the method (Alexander 1997), where data variance is esti- 
mated directly from internal scatter within the bin and not from 
external error estimates. This method was also used by Kaspi et 
al. (2007) in their analysis of the HET reverberation results. 
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The original f77 code by Dr. Alexander was reimplemented 
in f90 and f2py to form part of the python time series pack- 
age (Schwarzenberg-Czerny, 2012). Similarly, our fortran sim- 
ulation code was called via python wrapper. In doing so we took 
care to preserve continuity of our random series between sepa- 
rate fortran calls. 

For each of 1000 simulations of light curves calculated 
for a fixed parameter set we analyzed the DCF obtained from 
the ZDCF algorithm. To decrease scatter of the DCF it was 
smoothed with the 3-point running median filter For each 
smooth DCF we determined its maximum for positive lags /l,„ 
and attempted to fit it with a gaussian centroid function of 4 pa- 



2pl 



The fit 



rameters pi,i = 1, • • ■ ,4 f{A) - po + p\ exp 

performed with the python scipy . optimize . curve_£it suc- 
ceeded on 98 % occasions. So although the centroid could have 
been calculated for the majority of the cases, we decided to use 
the maximum since it works always and in general did not bring 
problems. 

The results of simulations yielded histograms of the max- 
imum lag Am and the centroid lag Ac s p4 statistics (Fig. [TJ. 
From them we conclude that the maximum of DCF statistics A„ 
yields the more precise and robust estimate of the true delay. As 
far as the limited scope of our simulations permits to establish, 
the scatter of the A„, statistics around its peak is less than that of 
Ac- Additionally, within precision of our simulations we were not 
able to find any bias of the estimated delay (lag of the histogram 
peak). 

4. Discussion 

We consider in detail the prospects of the determining the dark 
energy distribution from the reverberation studies of high red- 
shift quasars, as proposed by Watson et al. (2011). Our moti- 
vation comes from the observational campaign with this aim, 
done with the SALT telescope. However, measuring the delay 
between UV lines and the continuum in high redshift quasars 
is a telescope time consuming project and so a careful analysis 
performed beforehand can help to optimize it. 

Our Monte Carlo simulations of the planned monitoring pro- 
gram of high redshift AGN show that the success of the moni- 
toring crucially depends on three factors: 

- the choice of the mathematical method to determine the de- 
lay 

- the measurement accuracy 

- the choice of the emission line 

but it is also affected to some extend by other factors like the 
campaign duration and the distribution of the times of spectro- 
scopic and photometric measurements. These in fact are related 
to some extent to the choice of the quasars for monitoring. The 
simulated observational setup mimicked well what we can have 
for a single quasar with a large telescope, including sparse un- 
evenly distributed spectroscopic observations with gaps for the 
invisibility of the source. 

The choice of the method used to determine the delay does 
not have to be made in advance, but the estimate of the accuracy 
achievable in the whole project depends directly on this decision. 

We have checked three methods: ICCF (Interpolated Cross 
Corelation Function), ZDCF (Z-transformed Discreet Cross 
Correlation Function) and fitting. The first method almost 
always gave the largest error The ZDCF was much better, al- 
though occasionally it also failed due to development of the 
strong low values tail of the time delay histogram. The results 



were somewhat better, if a gaussian centroid fit was used in- 
stead of a maximum of the histogram, but in a few per cent of 
cases even the centroid method failed. The method was usu- 
ally the most stable one, although it showed a tail at large delays. 
Therefore, it seems that the best strategy for a real observation is 
to use both ZDCF andx^, and compare the results. 

The accuracy of the photometric and spectroscopic measure- 
ments is particularly important if the variability of the quasar is 
low (observed F„a,- less than 0.2 in a five year campaign). If the 
campaign last much longer and/or variability is much stronger 
then the measurement errors are less important, since the red 
noise variability in the lightcurves will dominate the measure- 
ment scatter. 

The choice of the emission line for monitoring is absolutely 
crucial. Our simulated results given in Tables [T]and |2]were ob- 
tained assuming that there is no strong intrinsic variability in the 
line itself, and the observed variability is due to reprocessing 
of the variable continuum, broadened with a gaussian response 
of the BLR. Indeed, it seems that Low Ionization Lines do not 
show intrinsic variability, and this was the source of the suc- 
cess of the Hy6 monitoring. We plan to use another LIL - Mg II 
2978 A which is also likely just to reflect the nuclear variabil- 
ity. However, if one wants to use the CIV line, which is actually 
necessary for monitoring quasars at redshifts z > 2, a problem 
appears. CIV lines belongs to the HIL class, and so comes from 
a different part of the BLR. As evidenced by the results of the 
HET campaign (Kaspi et al. 2007), in most cases Fuar of the 
CIV line was much higher than F^,,- of the continuum, implying 
a strong intrinsic variability of the line, and only one tentative 
measurement was obtained. 

4.1. Mgll vs. CIV line properties and prospects of monitoring 

Mg II line properties correlate well with the properties of the 
Balmer lines in large quasar samples (Salviander et al. 2007; 
McGin et al. 2008; Shen et al. 2008; Wang et al. 2009; Shen & 
Liu 2012). The line was succesfully used for black hole mass de- 
termination from a single epoch method by a number of authors 
(McLure & Jarvis 2002; Vestergaard & Osmer 2009, Shen et al. 

2011) . 

Woo (2008) showed fliat for Mg II line the variability level 
is in general comparable to the continuum variability. However, 
one out of 5 studies quasars (the object CTS300) actually 
showed much higher variability in Mg II line (rms of 17 % ) than 
in 3000 A continuum (rms of 7%). The reason for this variabil- 
ity is unclear since otherwise the source have similar parameters 
to the four other objects. 

CIV line properties poorly correlate with Balmer lines in a 
given object (Baskin & Laor 2005; Netzer et al. 2007 Shen & Liu 

2012) although some initial studies did not indicate considerable 
problems (Vestergaard 2002, Vestergaard & Peterson 2006). We 
discuss this issue in detail below. 



4.2. The problem of HET monitoring 

In the reverberation project for high redshift quasars with the 
HET telescope (Kaspi et al. 2007) the obtained results were in- 
conclusive, with only one measurement of a tentative delay of 
595 days in S5 0836-1-71. Here we reproduce the key properties 
of the monitoring: (i) objects were observable for six months 
only; typically (not always) 3 measurements were done, sepa- 
rated by 2.5 months, (ii) for S5 0836-1-71 the continuum vari- 
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Fig. 3. The predicted time delays in 1000 realizations of the 
model for the use of method, with the observational setup 
for HEX monitoring and assumed 660 day delay inserted in sim- 
ulations. There is a considerable pile in the last two bins of the 
delay histogram at the longest time delay which increases with 
increasing intrinsic variability from 5 % to 26%. 

ability amplitude was Fyai = 0.124, while the mean uncertainty 
for continuum measurements was 2%, and 10% for the CIV line. 

We calculate the predicted delay as a function of the CIV 
variability level, considering that the line shows also intrinsic 
variability in addition to the response to variable continuum. 

We model the intrinsic line variability by formally increas- 
ing its measurement error, and we performed computations for 
such setup using the method, which usually give the best 
results. The variability level of the continuum is only slightly 
lower than in our model A, low variability case (see Table [1]), 
but the spectroscopic and photometric covering is much lower. 
For the F„„(CIV) = 0.233 the recovered delay is 805 + 265. 
One sigma error is large (see Fig.O, and the upper limit is not 
well constrained, with 23.7 % of the delays stored in the last 
two histogram bins. The main problem here is not the spectro- 
scopic coverage: without the large intrinsic variability we would 
recover the delay 650 + 36 days for the same observational pat- 
tern. 

This means that CIV in general would require unrealisti- 
cally long observations for the intrinsic red noise variability of 
the continuum to dominate over the intrinsic variations of CIV. 
Therefore, a better way would be still to use the Mg II line and 
do the line observations in the near-IR band. 

4.3. Consequences for cosmology 

Our simulations illustrate the feasibility of testing the dark en- 
ergy distribution using quasar reverberation method. We plan to 
determine the observed monochromatic luminosity directly from 
photometry, and the intrinsic luminosity in our method is ob- 
tained from the relation 

log RBLRlMgll] = 1 .495 + 0.5 logL44,2897, (4) 



in the quasar rest frame. Here we recalculated the relation from 
Czerny & Hryniewicz (20 11) to the Mg II line wavelength, this 
relation is based on assumption of the disk effective temperature 
Teff = 1000 K at the onset of the BLR, due to the dust formation 
in the disk atmosphere. If the dust properties are indeed independ 
of redshift, there is no error in the constants in this formula. 

The size of the BLR is measured from the line-continuum 
time delay, and the accuracy is given in Tables[T]and|2] Adopting 
a realistic setup (model A, method) we thus gave 5% ac- 
curacy in photometry and 17% accuracy in time delay, which 
translates to the 17% error in the luminosity distance. If we as- 
sume that much better accuracy can be actually achieved (model 
B, x^ method) the final error reduces to 3%. 

It is more convenient to express it using the distance modulus 
(in magnitudes) defined as 

fi = m-M = 5logdL(z) + 25 , (5) 

where diiz) is the luminosity distance in megaparsecs and m 
and M are respectively the apparent and absolute luminosity of 
the source. The accuracy achievable in our method for a single 
source translates to 0.34 mag in the first case and 0.06 mag in the 
second case, for a single objects. By monitoring more objects we 
can increase these accuracy. 

In order to see whether this accuracy will likely bring inter- 
esting results we compare it with several cosmological models. 
For a given cosmological model one can compute the distance 
modulus as a function of redshift and compare it with observa- 
tions. In Figure m we plot residua of distance moduli with re- 
spect to the fiducial ACDM model, /imod - A'fid, assuming null 
spatial curvature (generalization to non-flat models is straight- 
forward). We show two wCDM models with a constant equa- 
tion of state of dark energy w = -0.8 and w - -1.2, and one 
with variable w(z) - wo + Waz/{1 + z) (the 'itioWaCDM' model). 
The fiducial cosmological parameters used were the central val- 
ues obtained by Mehta et al. (2012): = 0.28 and Hq = 70 
km s 'Mpc ' in case of the ACDM and wCDM models, and 
Qm = 0.272, Ho = 71km s-'Mpc"', wo = -1.02, iy„ = -0.26 for 
the WQWaCDM model. 

Since the expected differences between different cosmologi- 
cal models are of the order of 0.1 mag, it is possible to test the 
cosmological models with quasars at z ~ 1, if the measurement 
accuracy and/or the number of sources is high. 
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